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The El strength functions and radiative capture cross sections for several compound Sn isotopes, including unstable 
'^^Sn and '^"Sn, have been calculated using the self-consistent microscopic theory. In addition to the standard RPA or QRPA 
approaches, the method includes the quasiparticle-phonon coupling and the single-particle continuum. The results obtained 
show that the phonon contribution significantly affects the pygmy dipole resonance, which is of particular relevance for 
a proper description of the radiative neutron capture. The impact of the phonon coupling on the pygmy dipole resonance 
and the radiative neutron capture cross sections increases with the (A^ - Z) difference. For example, in the (0-10) MeV 
interval the full theory gives 17% of the energy-weighted sum rule for '^"^Sn and 2.8% for '^"'Sn, whereas the continuum 
QRPA approach gives 5.1% and 1.7%, respectively. These results indicate the importance the self-consistent calculation 
can have, especially when applied to neutron-rich nuclei of astrophysical interest. The comparison with the widely-used 
phenomenological Generalized Lorentzian approach shows that the (Q)RPA approach gives an increase in the neutron capture 
cross section by a factor of 2 for '^^Sn and a factor of 10 for '^"Sn and that the inclusion of the phonon coupling still increases 
these cross sections even furhter, by a factor of 2-3. 
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1. INTRODUCTION 

Inclusion of the coupling of single-particle degrees of 
freedom with phonon degrees (in short, phonon coupling 
or PC), in addition to the standard Random Phase Approx- 
imation (RPA) for magic nuclei or the Quasi-particle RPA 
(QRPA) for non-magic nuclei, is a main line of the recent 
development of microscopic nuclear theory. In particular, 
the Extended Theory of Finite Fermi Systems (ETFFS) (|T] 
includes both PC and the single-particle continuum, the 
latter being absolutely necessary for nuclei with the nu- 
cleon separation energy near zero. This approach has been 
recently generalized to include pairing in the quasiparticle 
time-blocking approximation ||2l, (ETFFS(QTBA)). This 
and other approaches developed with the PC have been 
supplemented with consideration of the self-consistency 
between the mean field and the effective interaction |I3]|4] 
|5|. On this basis, it was possible to perform a transition 
from the two sets of parameters used (the first one for the 
effective force and the second one for the mean field) to one 
unique set based on Skyrme force. Both these improvements 
-the single-particle continuum and the self-consistency- are 
of great importance, for astrophysics as well as for nuclear 
data evaluations, since they provide a reliable framework to 



calculate the structure of exotic nuclei, especially those with 
great neutron excess and/or with the nucleon separation 
energy close to zero. 

The PC role in the description of giant resonances in 
stable nuclei is well known: it explains approximately 50% 
of the observed widths, their gross structure and sometimes 
even the fine structure 1 1 1. However, the direct influence of 
the PC on giant resonances in unstable nuclei is less studied. 
It is clear that the PC role should be important for the widths 
here as well. But so far, the characteristics of the giant 
resonances in unstable nuclei have only been systematically 
studied within the (Q)RPA |6,7 8| approach. 

The PC plays an important role in the determination 
of the so-called pygmy dipole resonance (PDR), which is 
the low-energy part of the giant dipole resonance (GDR) 
that can exhaust about 1 - 2% of the energy-weighted sum 
rule (EWSR). At present there is no clear understanding 
around some important questions related to PDRs, see EH] . 
This resonance is also known to be of particular relevance 
in the description of the radiative neutron capture ll8][T0l. 
However, the PDR has essentially been treated on the basis 
of phenomenological models up to now. The question arises 
for exotic nuclei where the phenomenological approach 
may fail because of the specific features of the PDR in such 
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nuclei ||6llHl- The answer is the same as for the problem of 
giant resonances in unstable nuclei and is quite obvious: it 
is necessary to use a reliable theory for such nuclei, i. e., 
as discussed above, a self-consistent theory which accounts 
for the PC and the single-particle continuum in addition to 

the standard (Q)RPA. 

In our previous works |4, 9, 11] we developed the self- 
consistent version of the ETFFS(QTBA) using a discretiza- 
tion procedure for the single-particle continuum with dif- 
ferent kinds of Skyrme forces, including SLy4, where the 
velocity-dependent parts of forces were considered in the 
local approximation (we called it DTBA). The consequence 
of such a simplification is the necessity to renormalize the 
interaction in order to obtain the spurious state at zero 
energy. 

The main goals of the present work are: i) to study the 
PDR and GDR consistently in stable '^'^Sn and unstable 
'"'^Sn, '^"Sn isotopes using the variant (DTBA) of the 
microscopic self-consistent version of the ETFFS; ii) to 
investigate the impact of the PC on the radiative neutron 
capture cross section. 

2. SELF-CONSISTENT CALCULATIONS OF THE 
PDR AND GDR 

We use the SLy4 Skyrme force JTS) which proves to be 
rather successful in description of ground states and some 
excited states within the (Q)RPA. The ground states are 
calculated within the Hartree-Fock-Bogolyubov approach 
using the HFBRAD code |13|. The residual interaction 
for the (Q)RPA and the following QTBA calculations is 
derived as the second derivative of the Skyrme functional 
(see details in |4|). 

Fig. 1 shows the photoabsorption cross sections for 
'^"^Sn and '^"Sn nuclei. A more detailed discussion on 
'^^Sn can be found in [llj. One can see that the phonon 
contribution is noticeable for the PDR and it is increased 
with the (N-Z) difference growth. For example, in the (0-10) 
MeV interval the full theory gives 17% of EWSR for '^"Sn 
and 2.8% for '^"*Sn, whereas within the QRPA approach we 
have 5.1% and 1.7%, respectively. 

Fig. 2 shows the comparison of the El strength func- 
tions obtained in the framework of our self-consistent 
DTBA version with the (Q)RPA calculations. We used a 
normalization procedure here. Namely, the DTBA strength 
has been folded with a Lorentzian of 1 MeV width for 
DTBA and 2.7 MeV for RPA in order to broaden the GDR 
to a typical width of 4 MeV, as observed experimentally 
for the stable nuclei Although the DTBA peak energy of 
^^^Sn is predicted at a relatively low energy (13.5 MeV) 
with respect to experimental data (16.1 MeV), the PDR is 
correctly reproduced around 9.8 MeV [14|. 

In Fig. 2, our results are also compared with those ob- 
tained with the phenomenological Generalized Lorentzian 
strength function |15|. The Lorentzian approach is widely 
used for practical application, though it suffers from short- 
comings of various sorts. On the one hand, the location of 
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Fig. 1. Photoabsorption cross sections for '^^Sn (top) and ''"Sn 
(below) calculated within the self-consistent ETFFS (DTBA) 
without (dotted curves) and with (solid curves) the PC. 



the GDR maximum energy and width remain to be predicted 
from some underlying model for each nucleus. For many 
applications, these properties have often been obtained from 
a droplet-type of model or from experimental systematics 
1 10|. In addition, the Lorentzian model is unable to predict 
the El strength at energies below the neutron separation 
energy. Different parametrizations or functional forms (in- 
cluding in particular an energy- and temperature-dependent 
width) have been proposed (see e.g. |15|) to reconcile 
experimental data in the photon or radiative neutron capture 
channels, but none of the proposed closed forms can nowa- 
days explain the various trends observed at low energies. 
Most of all the Lorentzian approach cannot provide any 
predictions on the low-energy PDR, neither on its presence, 
nor on its characteristics. For this reason, it is of particular 
interest to analyze to what extent our predictions based on 
self-consistent microscopic models differ from those used 
in practical applications. One can see that the difference 
between [15|, on the one hand, and both the QRPA and 
DTBA, on the other hand, increases with the (A^ - Z) 
difference. As our preliminary calculations have shown, 
the main reason is that the A-dependence of the integral 
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characteristics of neutron-rich nuclei differs strongly from 
the usual ones. This fact demonstrates the non-applicability 
of the standard systematics (e.g. £0 - SOA"'''^) for unstable 
nuclei. 



3. CALCULATIONS OF RADIATIVE NEUTRON 
CAPTURE CROSS SECTIONS 
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Fig. 2. Strength functions for the isovector 1" states in '-''Sn, 
'^^Sn and ""Sn nuclei. The solid, dashed and dotted curves are 
for the DTBA, (Q)RPA and Kopecky-Uhl I10II15II approaches, 
respectively 
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Fig. 3. Radiative neutron capture cross sections obtained with the 
y-ray strength functions of '^^Sn, '^^Sn, '^"Sn shown in Fig. 2 



The radiative neutron capture cross sections have been 
calculated using the reaction code TALYS [16|. The DTBA 
strength functions with and without PC have been included 
in the calculation of the electromagnetic de-excitation trans- 
mission coefficients. The final radiative neutron capture 
cross sections calculated with the strength functions of 
Fig. 2 are shown in Fig. 3. We show the cross sections 
leading to the '^'^Sn, '^^Sn and '^"Sn compound nuclei, 
respectively. Note that such reactions are of interest for 
the r-process nucleosynthesis, and illustrate the impact that 
the low-lying strength can have on cross sections, and 
consequently on reaction rates. As seen in Fig. 3, including 
the PC increases the cross section by a factor of 2-3. In 
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all cases the trend is clear: the cross sections follow on the 
whole the strength functions shown in Fig. 2, and the DTBA 
with an extra low-lying strength gives a final reaction cross 
section (and Maxwellian-averaged rate around T - 1 - 2 ■ 
\{fK) about 3 times larger compared with the prediction 
obtained with the HFB+QRPA calculation |7|. For the case 
of the compound '^"^Sn, we see in Fig. 2 that the difference 
between the DTBA and QRPA strength functions is much 
smaller than for '^^Sn and '^"Sn. It is necessary also to note 
here that the theoretical consideration for '^"^Sn is based in 
a sense on the known experiment, this is especially true for 
the Kopecky-Uhl approach (see below). For these reasons, 
the appropriate cross sections in Fig. 3 are insensitive by 
sight to the three strength functions. 

In Fig. 3 we also compare the cross sections with the 
widely used Generalized Lorentzian strength function IITOl 
[TSll . We have found that for the stable '^"^Sn nucleus the 
cross section is almost identical to the one obtained with 
the ETFFS approach, although the strength functions can 
differ significantly at energies below the neutron separation 
energy. For neutron-rich isotopes, the QRPA strength gives 
a significant increase in the cross sections with respect to 
those obtained with the Lorentzian function of LlSi . namely 
by a factor of 2 for '^^Sn and a factor of 10 for '^°Sn. 
This confirms the results of |7|. Moreover, the inclusion 
of PC increases the cross sections for the unstable '^^Sn 
and ^^"Sn by a factor of 2-3 at the energies around 100 
keV, as illustrated in Fig. 3. These results demonstrate the 
necessity of self-consistent approaches for the calculation 
of the radiative capture cross sections of exotic neutron-rich 
nuclei. 



4. CONCLUSION 

The electric dipole strength function has been calculated 
on the basis of the ETFF(DTBA) model which simul- 
taneously takes into account the RPA or QRPA config- 
urations, the more complex lplh(8iphonon configurations 
and the single-particle continuum. These are the necessary 
ingredients for a consistent study of the giant and pygmy 
resonances. 

It was demonstrated that the DTBA approach predicts a 
relatively large low-lying strength compared to the (Q)RPA 
one around 10 MeV in the neutron-rich Sn isotopes. The 
radiative neutron capture cross sections were calculated 
with the DTBA and (Q)RPA strength functions and shown 
to be sensitive to the predicted low-lying strength. In par- 
ticular, the calculations demonstrated that including the PC 
leads to a significant increase in the reaction cross section. 
The comparison with the Generalized Lorentzian approach 
confirmed the non-applicability of the phenomenological 
approaches for neutron-rich nuclei. 

The ETFFS(DTBA) approach is believed to provide a 
more complete and coherent description of the gamma-ray 
strength function than the previous models used so far For 
astrophysics applications in particular, such calculations are 
highly recommended for a more reliable estimate of the 
electromagnetic properties of exotic nuclei. 
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